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An exact analytical solution of the statistical multifragmentation model is found in thermody- 
namic limit. The system of nuclear fragments exhibits a 1-st order liquid-gas phase transition. The 
peculiar thermodynamic properties of the model near the boundary between the mixed phase and 
the pure gaseous phase are studied. 
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Nuclear multifragmentation is one of the most inter- 
esting and widely discussed phenomena in intermediate 
energy nuclear reactions. The statistical multifragmenta- 
tion model (SMM) (see and references therein) was 
recently applied to study the relationship of this phe- 
nomenon to the liquid-gas phase transition in nuclear 
matter [^-^. Numerical calculations within the canon- 
ical ensemble exhibited many intriguing peculiarities of 
the finite multifragment systems. However, the inves- 
tigation of the system's behavior in the thermodynamic 
limit is still missing. Therefore, there is no rigorous proof 
of the phase transition existence, and the phase diagram 
structure of the SMM is unknown yet. In the present 
paper an exact analytical solution of the SMM is found 
within the grand canonical ensemble. The self-consistent 
treatment of the excluded volume effects is an important 
part of our study. 

The system states in the SMM are specified by the 
multiplicity sets {uk} {n-k = 0, 1, 2, ...) of fc-nucleon frag- 
ments. The partition function of a single fragment with 
k nucleons is 0: oJk = V Zk , where 

k = 1,2,..., A {A is the total number of nucleons in the 
system), V and T are, respectively, the volume and the 
temperature of the system, m is the nucleon mass. The 
first two factors in ujk originate from the non-relativistic 
thermal motion and the last factor, Zk, represents the 
intrinsic partition function of the /c-fragment. For fc = 1 
(nucleon) we take zi = 4 (4 internal spin-isospin states) 
and for fragments with fc > 1 we use the expres- 
sion motivated by the liquid drop model (see details in 
Ref. 0): Zk = exp(— /fc/T), with fragment free energy 
fk = -[Vo + T^/€o]k + a(T)k^^^. Here Wo = 16 MeV is 
the bulk binding energy per nucleon, T^/eq is the con- 
tribution of the excited states taken in the Fermi-gas 
approximation {cq = 16 MeV) and a{T) is the surface 



tension which is parameterized in the following form: 
a{T) ^ ao[{T^ - T^)/{T^ + T^)f'^, with - 18 MeV 
and Tc = 18 MeV (cr = at T > T^). The canonical par- 
tition function (CPF) of nuclear fragments in the SMM 
has the following form: 



(1) 



In Eq. (^) the nuclear fragments are treated as point- 
like objects. However, these fragments have non-zero 
proper volumes and they should not overlap in the co- 
ordinate space. In the excluded volume (Van der Waals) 
approximation this is achieved by substituting the to- 
tal volume V in Eq. (|l|) by the free (available) volume 
Vf = V — bY^^kuk, where b = 1/po (po = 0.16 fm~^ 
is the normal nuclear density). Therefore, the corrected 
CPF becomes: Za{V,T) = Z'^^iV - bA,T). The SMM 
defined by Eq. (0) was studied numerically in Refs. P-p|. 
This is a simplified version of the SMM, e.g. the symme- 
try and Coulomb contributions are neglected. However, 
its investigation appears to be of principal importance 
for studies of the liquid-gas phase transition. 

The calculation of Za{V,T) is difficult because of the 
constraint kuk — A. This difficulty can be partly 
avoided by calculating the grand canonical partition 
function: 

oo 

Z{V,T,p) = Y,e^v(tif^ZA{V,T)Q{V ~bA) , (2) 

where chemical potential /i is introduced. The calcu- 
lation of Z is still rather difficult. The summation over 
{n/c} sets in Za cannot be performed analytically because 
of additional A-dependence in the free volume V} and the 
restriction Vf > Q. The problem can be solved M by in- 
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troducing the so-called isobaric partition function (IPF) 
which is calculated in a straightforward way: 



i'(s,T, ^) 



dV exp(-sV) Z{V,T,^l) 



J'{s,T,fi) 



(3) 



where 



T{s,T, fi) 



3/2 



+ f;fc3/2exp 

fe=2 



zi exp 



/i — s6T 
T 



(iy - sbT)k ~ Gk^l^ 



T 



(4) 



with V = ^ + Wo + T^/eo. In the thermodynamic limit 
V ^ oo the pressure of the system is defined by the 
farthest-right singularity, s*{T,fj,), of the IPF Z{s,T,^) 
(see Ref. M for details): 



T lim 



In Z(y,T,^) 



Ts*(r,M). (5) 



The study of the system's behavior in the thermodynamic 
limit is therefore reduced to the investigation of the sin- 
gularities of Z. 




FIG. 1. Phase diagram in the (T, p)-plane. The mixed 
phase and pure gaseous phase boundary is shown by the solid 
line. The pure liquid phase (shown by crosses) corresponds 
to the fixed density p = pa- Point C is the critical point, at 
T > Tc only the pure gaseous phase exists. 

The IPF (|) has two types of singularities: 
1) the simple pole singularity defined by the equation 
Sg(T, /i) — T{sg,T, /i) ; 2) the singularity of the function 
T itself at the point s/(T, /i) = v/Th where the coefficient 
in linear over k terms of the exponent in Eq. (^ is equal 
to zero. 

The simple pole singularity corresponds to the gaseous 
phase where pressure pg = Tsg is determined by the tran- 
scendental equation: pg{T,fj,) — TT{pg/T,T, jj,). The 



singularity s/ (T, fj.) of the function !F defines the liquid 
pressure: pi{T,n) = Tsi{T,Ijl) = v/b. Here the liquid is 
represented by an infinite fragment with k ^ oo. 

In the region of the (T, /x)-plane where ly < bpg (T, fi) 
the gaseous phase is realized (pg > pi), while the liquid 
phase dominates at v > bpg (T, /i) . The liquid-gas phase 
transition occurs when the two singularities coincide, i.e. 
Sg{T,fi) — si(T,fj,). As J- in Eq. (^ is a monotonously 
decreasing function of s the necessary condition for the 
phase transition is that the function J-' is finite in its sin- 
gular point si. This condition requires (t{T) > and, 
therefore, T < Tc- Otherwise, J-{si,T,ii) — oo and the 
system is always in the gaseous phase as Sg > si. 

At T < Tc the system undergoes a 1-st order phase 
transition across the line /i* = IJ*{T) defined by the con- 
dition of coinciding singularities: si = Sg . The baryonic 
density p in the liquid and gas phases is given by the 
following formulae, respectively: pi = {dpi/dpt)rp = 1/6, 
Pg = {dpg/dp)rj, = / ( 1 + ^^id ) , whcrc the function p^d 
is the density of point-like nuclear fragments with shifted, 
p, p — bpg, chemical potential: 
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A similar expression for pg within the excluded volume 
model for the pure nucleon gas was obtained in Ref. 
The phase transition line p*{T) in the (T, /i)-plane cor- 
responds to the mixed liquid and gas states. This line 
is transformed into the finite mixed-phase region in the 
(T, p)-plane shown in Fig. 1. The baryonic density in 
the mixed phase is a superposition of the liquid and 
gas baryonic densities: p = Xpi -|- (1 — X)pg , where A 
(0 < A < 1) is a fraction of the system's volume occu- 
pied by the liquid inside the mixed phase. Similar lin- 
ear combinations are also valid for the entropy density s 
and the energy density e with {i = l,g) Si = (dpi/dT)^, 
e, = T (dpi/dT)^ + p idpjdp)j, - p,. 

Inside the mixed phase at constant density p the pa- 
rameter A has a specific temperature dependence shown 
in Fig. 2: from an approximately constant value p/ Po 
at small T the function A(T) drops to zero in a nar- 
row vicinity of the boundary separating the mixed phase 
and the pure gaseous phase. Such an abrupt decrease 
of A(r) near this boundary causes a strong increase of 
the energy density as a function of temperature. This is 
evident from Fig. 3 which shows the caloric curves at dif- 
ferent baryonic densities. One can clearly see a plateau- 
like behavior. As a consequence this leads to a sharp 
peak in the specific heat per nucleon at constant density. 
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Cp(T) = {d£/dT)p/p , shown in Fig. 4. 




P/Po= 1'3 



FIG. 2. Volume fraction A(T) of the liquid inside the 
mixed phase is shown as a function of temperature for fixed 
nucleon densities p/po ~ 1/6, 1/3, 1/2, 2/3, 5/6 (from bottom 
to top). 



The peculiar thermodynamic property of the consid- 
ered model is a finite discontinuity of Cp(T) at the bound- 
ary of the mixed and gaseous phases seen in Fig. 4. This 
finite discontinuity is due to the fact that X{T) = 0, but 
{dX/dT)p ^ at this boundary inside the mixed phase 
(see Fig. 2). It should be emphasized that the energy 
density is continuous at the boundary of the mixed phase 
and the gaseous phase, hence the sharpness of the peak 
in Cp is entirely due to the strong temperature depen- 
dence of A(T) near this boundary. Therefore, in contrast 
to the expectation in Refs. ||,p|, the maximum value of 



Cp remains finite and the peak width in Cp(T) is not zero 
in the thermodynamic limit considered in our study. 




T (MeV) 

FIG. 4. Specific heat per nucleon as a function of temper- 
ature at fixed nucleon density p/po = 1/3. The dashed line 
shows the finite discontinuity of Cp{T) at the boundary of the 
mixed and gaseous phases. 

In conclusion, the simplified version of the SMM is 
solved analytically in the grand canonical ensemble. The 
progress is achieved by reducing the description of phase 
transitions to the investigation of the isobaric parti- 
tion function singularities. The model clearly demon- 
strates a 1-st order phase transition of the liquid-gas type. 
The considered system has peculiar properties near the 
boundary of the mixed and gaseous phases seen in the 
behavior of the specific heat. 
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FIG. 3. Temperature as a function of energy density per 
nucleon (caloric curve) is shown for fixed nucleon densities 
p/Po = 1/6, 1/3, 1/2, 2/3. 
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